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Abstract 

In this paper, we develop a new integral representation for the solu- 
tion of the time harmonic Maxwell equations in media with piecewise 
constant dielectric permittivity and magnetic permeability in . This 
representation leads to a coupled system of Fredholm integral equations 
of the second kind for four scalar densities supported on the material 
interface. Like the classical Miiller equation, it has no spurious res- 
onances. Unlike the classical approach, however, the representation 
does not suffer from low frequency breakdown. We illustrate the per- 
formance of the method with numerical examples. 



1 Introduction 

In our previous paper we introduced a new representation for the time 
harmonic Maxwell equations in M'^, based on two scalar densities defined 
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on the surface bD of a smooth bounded region D. This bounded region 
may have several components, but we assume that its complement is 
connected. We refer to these densities as generalized Debye sources, since 
they generalize to arbitrary geometries the classical formalism of Lorenz, 
Debye and Mie that is limited to the sphere. We also showed in [1] that 
the problem of scattering from a perfect conductor can be reduced to the 
solution of a coupled pair of Fredholm boundary integral equations of the 
second kind. This system of equations was shown to be invertible for all 
non-zero wave numbers in the closed upper half plane. Moreover, in the 
case that all components of bD are simply connected, we showed that this 
system of equations does not suffer from a phenomenon called "low frequency 
breakdown" . 

Here, we develop an integral equation for the case of dielectric (interface) 
boundary conditions and extend the analysis of low frequency breakdown to 
the multiply connected case. We use, almost exclusively, the representation 
of fields in the language of forms, see We let 

3(x, t) = i{x)e-'^^ N{x, t) = v{x)e-'^\ (1) 

where H is a 1-form representing the electric field and AT is a 2-form repre- 
senting the magnetic field, that is -E -f-T- ^ and H ^ r/. Faraday's law and 
Ampere's law (the curl equations) take the form: 

d$, = iujfiT] d*r} = —iuje$,. (2) 

For Lo ^ 0, these equations imply the divergence equations, which take the 
form: 

d*$, = d'n = 0. (3) 

Together, ([2]) and ^ give an elliptic system for the pair {$,,ri). In the di- 
electric/interface case, our representation for the THME(A;) applies within 
the bounded components of D as well as in the exterior domain. The per- 
mittivity e and permeability fj, inside each component may be distinct from 
the corresponding values in the exterior. 

More precisely, following the discussion in [7], we assume that in each 
component of \ bD, the EM-parameters are piecewise constant with 

_ .(T _ .a' . ^. 

e = e + i— fi = fi + i—. (4) 

UJ u 

Here, a and a' are non-negative numbers, while ?and Jl are positive numbers. 
The complex numbers we and ufi lie in the closed upper half plane. We, 
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therefore, assume that the arguments of their square roots he in the interval 
[0, §]. Hence if 

k = uj^feii = ^/uJe^/ujJi, (5) 

then we can assume that 

< arg/c < vr provided < argw < vr. (6) 

We restrict our attention here to the standard case where u G M+. We call 
the system of equations ([2]) and ([3]) the THME(/c). 

The EM-parameters of = \ D, are denoted (ei,^i), with ki = 
u)^e\[i\. If the components of D are {-Di, . . . , -Da?}, then the EM-parameters 
for Dj are (eoj) /^Oj)) with k^j = Lv^eojfiQj. The total EM-field in Q is written 
in the form 

Here is an arbitrary solution to the THME(A;i) defined in fi, and 

{^^'^,rj^^) is an outgoing solution to THME(/ci), selected to insure that 
(^*°*, ?7*°*) has the correct behavior on bD. The superscript is usually omit- 
ted in the sequel. In the perfect conductor case, the fields in the bounded 
components of D are zero, while in the dielectric case the total field is each 
bounded component is just a scattered field which satisfies the appropriate 
variant of Maxwell's equations. 

As noted, the scattered field, {$^^, V^'^), is assumed to satisfy the outgoing 
radiation condition in fi. For the electric field, this reads: 

^^d^ - ^k^ = O ^-^y ^ = ^(r)' 

where x = The same condition is satisfied by -k^r/, where +3 is the 

Hodge star operator acting on forms defined in M^. It is a classical result 
that if (^, Tj) solves the THME(A;) for k ^ 0, with non-negative imaginary 
part, and one component is outgoing, then so is the other. When uj = (or 
k = 0) the equations for ^ and r] decouple; the divergence equations, ([3]), are 
no longer a consequence of the curl equations, but are nonetheless assumed 
to hold. 

In this paper, we continue our study of the representation of solutions to 
the time harmonic Maxwell equations in terms of the scalar Debye source 
densities, which we denote by (r, g), supplemented, in the case that bD is 
of genus (7 > 0, by the 25(-dimensional space of harmonic 1-forms j fj. We 
first show how to use our representation to solve the time harmonic Maxwell 
equation when D is a dielectric with piecewise constant fi, e and a, under the 
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physical boundary condition that the tangential components of the ^ and ij 
fields are continuous across bD. As before we obtain a system of Fredholm 
equations of second kind, which does not suffer from either spurious, interior 
resonances or low frequency breakdown. The key to the good behavior as 
A; — )• 0, is that the scalar Debye sources (r, q) are most directly related to the 
normal components, along bD, of the $, and r] fields. Unlike the tangential 
components, which, at A; = 0, must satisfy a differential equation along bD, 
{dhDCt = 0) the normal components are not in any way constrained when 
/c = 0. Thus, as a means of parameterizing solutions of the THME(/!;), the 
normal components behave much better as A; — t- than does the tangential 
data. 

In the second part of the paper we give a detailed analysis of the low 
frequency behavior of our representation in the perfect conductor case. In 
the case that the genus of bD is non-zero, the dielectric problem is somewhat 
simpler than the perfect conductor. At A; = the solutions of the time 
harmonic Maxwell equations are harmonic fields: 

d*$^ = d^ = and d*r] = dr] = 0. (8) 

This means that these fields represent cohomology classes in H^^{Cl) and 
-ffjj^(il), respectively. Therefore, in the case of perfect conductors, the re- 
strictions of , and ^ar?"*" to bQ = bD each span subspaces of H^j^{bD) 
of half the total dimension. In the dielectric case, we need to consider the 
jumps of these fields across bD. For all wave numbers, including A: = 0, the 
harmonic projections of these jumps span all of H^j^{bD). Hence the topo- 
logical constraint that arises in the perfect conductor problem, complicating 
the low frequency behavior, is absent in the dielectric case. 

In Section [2] we review the generalities of the representation of solutions 
in terms of Debye sources, including the restrictions to the boundaries and 
the jump conditions. In Section [3] we show how to use this approach to 
represent solutions to the dielectric problem and prove the basic unique- 
ness results. In Section |4] we derive the boundary integral equations for 
the dielectric problem and show that they do not suffer from low frequency 
breakdown in either the simply connected, or non-simply connected cases. 
Finally, in Section [5l we present a new approach to the perfect conductor 
problem, when bD is not simply connected. Using this approach there is 
no low frequency breakdown as A; — )• 0, and the perfect conductor prob- 
lem for the THME(A;) gracefully decouples. Section [6] contains numerical 
experiments illustrating several of the results proved in earlier sections. 
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2 Debye Sources and Potentials 



We begin by reviewing the symmetric representation of solutions to the 
THME(A;) in terms of both potentials and anti-potentials. We assume, as 
discussed above, that the time dependence is e"*"^*, and that the permittivity 
e, the permeability fi, and the conductivity a are piecewise constant. We let 
k = ujy/ejl, which we take to have non-negative imaginary part. As in [5, 
we represent the solution to the THME(A;) by setting: 

= ^{ikcx- d<j)- (fO) r] = ^e{ike-d*-^ + dci), (9) 

where (/> is a scalar function, cx a one form, 6 a two form, and ^ = tpdV, a 
three form; a is the usual vector potential and (p the corresponding scalar 
potential, while 6 is the vector anti-potential and ijj the corresponding scalar 
anti-potential. There are many possible choices for scaling the coefficients of 
the various terms on the right hand side of ^ . An advantage of the scaling 
in ([9]) is that, in each sub-region, /x and e only appear as global multipliers, 
while the terms within the parentheses depend only on the wave number 
k = ijj^JJie. 

Assuming that all of the potentials solve the Helmholtz equation, A/3 + 
= 0, where A = — + dd*) denotes the (negative) Laplace operator 
in the correct form degree, for (^, r]) to satisfy the equations in : 

= iujiiT] d*r] = —iue^, (10) 

it suffices to check that (in the Lorenz gauge) 

d*a = -ik(p de = ik^. (11) 



2.1 Debye source representation 

We let gk{x,y) denote the fundamental solution for the scalar Helmholtz 
equation, with wave number k, which satisfies the Sommerfeld radiation 
condition: 

^ik\x-y\ 



9k{x,y) 



4tt\x — y\ 

For the moment we assume that D is connected, and let T = bD. All of the 
potentials can be expressed in terms of a pair of 1-forms j,m defined on 
r, which define electrical and magnetic currents. As T is embedded in 
these 1-forms can be expressed in terms of the ambient basis from M^, e.g., 

j = h{x)dxi + j2{x)dX2 + j3{x)dx3] (12) 
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we normalize with the requirement 



ir.j=j{n) = 0. (13) 
It is well known [6l [9] that the conditions (jlip are satisfied if 

a = y" gkix,y)[jiy) ■ dx]dA{y) (j) = ^ j gkix,y)dr *2 jiv) 

^ ' ^ (14) 

f dVx f 

= *3 9k{x,y)[m{y) ■ dx\dA{y) ^ = j gf,{x,y)dr *2 m{y). 

r r 

In the end, however, we do not use the currents j and m as the "funda- 
mental" parameters. In [3], we introduced the notion of generalized Debye 
sources, r, q, which are scalar functions defined on T : 

4"'^r *2 j = fdA 4"'^r *2 'ti = qdA. (15) 
ik ik 

From this definition, we see that rdA and qdA are exact and hence their 
mean values vanish on F, 



rdA = j qdA = (16) 
r r 

This is necessary for the conditions in (|15p to hold, and thus, for (^, rf) to 
satisfy the Maxwell equations. In terms of the generalized Debye sources: 

(k = j 9k{x,y)r{y)dA{y) ^ = dV^ J gk{x,y)q{y)dA{y). (17) 
r r 

Below, we derive equations for j and m in terms of these scalar potentials 
and, if needed, harmonic 1-forms. 

Provisionally, we let T^{e, fi,uj,T, j ,m) denote the fields defined by dH) 
and ()14p . Here + refers to the unbounded component of F"^ and — , the 
bounded component. In the sequel, j and m are usually taken to be func- 
tions of scalar sources, e.g. r and q, via the relations in (llSp . along with a 
possible harmonic components [j fj^mn)- The functional relationships be- 
tween the currents {j,m) and the Debye sources {r,q,jf^,mH) depend on 
the relationships between the currents themselves, which in turn depend on 
the particulars of the boundary value problem we are trying to solve. For 
example, in the perfect conductor problem, we take m = ★2J) and then 

j = ik[drRor - i^2drRoq] + 3hi (18) 
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where Rq is the inverse of the (negative) scalar surface Laplacian, 

^r,ou = {*2dr *2 dr + dr *2 dr*2)u, (19) 

restricted to functions of mean zero. Note that, in the hmit w — )• 0, only the 
harmonic components of the surface currents play a role. 

2.2 Mapping properties of the Debye source representation 

The regularity of the fields (^^ , r/^ ) is straightforward to describe in terms 
of the regularity of the Debye source data: {r,q,j^,TnH)- The harmonic 
components {jfj^mu) are always infinitely differentiable. If r and q belong 
to the L^-Sobolev space, H^{r), then the currents j and m, defined by (fTSj) 
belong to H'^'^'^iT). As we shall see, this remains true for the dielectric 
problem, even though the relationships amongst the currents and the scalar 
sources are somewhat different. 

For any k in the closed upper half plane, and real number s, the single 
layer potential defines bounded maps 

St : ^'(r) H'+hn) and 

^ ' , ^ (20) 

S- : H'iV) ^ H'+2(D) 

From these observations we conclude that, if {r,q) S H^{T), then 
a+,6'+ G H'+^n), ,0+,^+ G H'+'^n), and 

Taken together, these observations along with ^ show that if (r, q) G H^(T), 
then 

(,C+,r7+) G H'+^n) and G H'+^D). (22) 

In our applications the Debye sources (r, q) are determined by solving 
Fredholm equations of second kind, with kernels defined by elliptic pseudod- 
ifferential operators of order 0. Hence, if the data belong to H^{T), then so 
do the Debye sources, and therefore the fields in F'^ belong to H^~^2. This 
is precisely what one expects for the solution of an elliptic boundary value 
problem for a first order elliptic system. 

2.3 Boundary equations 

Following the convention in [3], we use *2^,t *2([*3^]t)) which correspond 
to n X E and n x H, respectively, to represent the tangential components. 
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and the scalar functions in^^ and Zn[*3?7] to represent the normal compo- 
nents, corresponding to n- E, and n - H. These limiting values are given by 
the integral operators: 



*2([*3»7^]t/Ve. 



1 /±m 
2 



+ 











q 
3 



'n 



(23) 



and 



«n[*3?7'^]/A/e 



1 f±r 
2\±q 



-Ko 








K3 



ikK2,n 



(r\ 
Q 
3 

\mj 



d I'M} 



(24) 



The operators Kq, Ki, K2^t, K2^n, and are defined in the appendix. 

In what follows, we allow the bounded domain D, and therefore its 
boundary, T to have several connected components, {Fi, . . . , Fat}. In all 
cases, we let M.r,o denote pairs of functions (r, q) defined on F so that 



rdA = J qdA = for I = 1,...,N. 
r, r. 



(25) 



Such functions are referred to as scalar Debye sources. We let ^^(F) denote 
the vector space of harmonic 1-forms on F, that is, the solutions of 



drcx. = dpo; = 0. 
Let pi be the genus of F/, and 

P = Pi-\ \-PN, 

the total genus of F. It is a classical theorem that 

dim-H^Ti) = 2pi. 

Thus dimWi(F) = 2p. 



(26) 
(27) 
(28) 
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3 Uniqueness for the Dielectric Problem 



We now apply our representation to the problem of several dielectric mate- 
rials separated by smooth bounded interfaces. We begin with the slightly 
simpler case of a single connected, bounded region. The bounded region is 
denoted by D; we assume that hD = F is connected and is the exterior 
region \ D. We let (eo,/io) denote the EM-parameters for D and (ei,/ii) 
denote the EM-parameters for O. Given a frequency a; 7^ from the closed 
upper half plane, we set 

ki = Uy/JUTi for / = 0, 1, (29) 

where < axgki < vr. We denote by (^^'"'^,''7^'"'^) the electromagnetic field 
in the exterior region corresponding to the exterior parameters and by 
{^ '^,T]^'^) the electromagnetic field in the interior region D corresponding 
to the interior parameters. 

The dielectric problem for (^, rj) involves the determination of these fields 
satisfying 

^P' - = jf and (z„r7+'^)t - {ir^rj-^ = - *2 ml^ (30) 

where (j'",m"^) are two 1-forms specified along V. Where we recall that n 
is the outward normal, relative to D along T, and that: 

*2 [inll) = {*3V)t- (31) 

Note that, in the case of a scattering problem, (j'°,m™) are the tangential 
components of the known incoming field. 

For I = 0, 1, we suppose that these solutions are defined by currents 
{jl,mi) defined on F. For Z = 0, 1 we let 

{$^'\V^'') = ^^{ei,fii,co,r,ji,mi) (32) 

denote the solutions to the THME(fc/) in \ F specified by these sources. 
Solutions defined by this ansatz automatically satisfy the outgoing radiation 
condition. 

We seek a solution of the dielectric problem given by 

^) = /•^~(^o,Mo,w,F,jo,mo) inD ^^^^ 
1 J'+(ei,//i,w,F,j\,mi) inn, 
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according to ([9]), (jl4p . Note, however, that we have four unknown 1-forms 
on r but only two 1-forms as data. We therefore, suppose that the inner 
and outer currents are related by a transformation of the form 

jo = \ —Uji and mo = . —Umi, (34) 
V eo V /"o 

where U is what we refer to as a clutching map. 

Definition 1. A clutching map is a linear isomorphism 

U :C^{bD;A^) ^C^{bD;A^), 

that satisfies the conditions: 

1. The map is complex symplectic: for a,l3 £ C^{bD; A^) we have 

j)aAp = ^UaAUp. (35) 

bD bD 

2. The map U preserves the harmonic forms ^(^{^{bD) = ^{^{bD). 

3. On the orthogonal complement of the harmonic 1-forms, 

n^bD)^ = dbDC\bD)®dljjC\bD;A'^), 
U reduces to the Hodge star-operator: 

^ \'H^{bD)^= *2- (36) 

The simplest example is to set: 

= *2j, 

For this choice, however, a mild form of low frequency breakdown occurs in 
certain non-generic cases. 

The harmonic 1-forms are invariant under the action of *2 • It is also the 
case that if a G Ti^{bD) and /3 = d^DU + d^j^vdA, then 



Q A /3 = y a A *2/3 = 0. (37) 

bD bD 

Prom these observations and condition 3, it is apparent that U is nothing 
more than a choice of hermitian symplectic isomorphism from 1-0- {bD) to 
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itself. If D has multiple components, {-Di, . . . j-Dat}, then this is done one 
component at a time, i.e. we choose maps lAi : H^ipDi) — t- H^ipDi), for 
/ = 1,... ,iV. 

With scalar Debye sources {ri,qi) and the requirement 

dr *2 Ji = ikiridA and drmi = ikiqidA, (38) 
conditions 2 and 3 on the map Vl imply that: 

Aiji = iuj{^/jr[eldrri - 60^^*2 dr?'o) (39) 

Aimi = iuj{y/meidrqi - MoW — *2 rfrQo)- (40) 

V /^i 

This means that we can define ji and mi much as in the perfect conductor 
case: 

ji = iuj{y/jr[e[drRori - eoW — *2 drRoro) +3h-, (41) 

V 

and 

mi = iuj{^fiieidrRoqi - fJ'0\ — *2 drRoqo) + mn- (42) 

V 

Here i?o is the partial inverse of the (negative) scalar Laplace operator on V 
and jfj,mH are the harmonic projections of Ji,mi, respectively. Clearly, 
Ji and mi are of order —1 in terms of the scalar sources. If w = 0, then j'l 
and mi are purely harmonic. 

In Theorem [T] below, we also make use of the jump relations 

^t^'' - ^7'' = *2 rni and [1^1]'^'% - {inr]~'% = V^di, (43) 



which follow for / = 0, 1 from equation (j23p . Following Miiller's argument, 
we now show that our parametrization satisfies a basic uniqueness require- 
ment for all non-zero frequencies in the closed upper half plane. 

Theorem 1. Assume that oj ^ 0, that ujfj,o,ujeo,u}fii,ujei have non-negative 
imaginary parts, and that Re(/io/eo) > 0. Let bD be connected and 

^ = {ro,qo,ri,qi,jfj,mH) 
be Debye source data defining solutions to THME(ki) in bD^, with 

Jo = \ —^Ji 0'''T'd 1^0 = \ —Umi, (44) 
V eo V /^o 
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for U a clutching map. If ij) given by (j33p satisfies: 

it - C = and (^r^^7+)^ - {ir,V')t = 0, (45) 
then the data, (?'0; 90) ^ij 91) j'/zi ^/f)) are also zero. 

Remark 1. For non-zero frequencies, the representation is unique for any 
choice of clutching map. After proving this theorem we consider what hap- 
pens if D has several components, and finally what happens when a; = 0. 

Proof. Since (^^'^, ^7^'^) is an outgoing solution to THME(A;i) it follows from 
Miiller's uniqueness theorem (Theorem 61 in[7j) that (^^'"^, ^7"*"'^) = (0,0) 
and (^"•°,T7-'°) = (0,0). To prove the theorem we need to show that 

,~ ~. { (^+'°,T7+'°) in J7 , , 

1 (r'\^-'^)mD. ^ ' 



also vanishes. The jump conditions (j43p imply that 

= -VJm)*2 mo, {inV^% = VeoJo) (47) 

and 

^7'^ = ^//iT*2 mi, {inr]'''^)t = -\/eiJi- (48) 
The boundary current relation (j44p then implies that 

^+'° = -Vm*2Z^mi, {ir.V^'% = V^Uj^. (49) 

Stokes theorem shows that: 



bD 



D bD 



(50) 



bD 



bD 



We use f^-^ to emphasize that this is the pairing of a 2-form with a 2-cycle, 
where bX is oriented as the boundary of X. If = n Bji{0), then a 
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second application of Stokes theorem gives 

Qn bn bBn 

bBR 



hn bBR 



(51) 



The orientation of T as the boundary of 17 is opposite to that as the boundary 
of D, and *2 is a point-wise isometry, hence ()35p shows that these boundary 
contributions are of equal magnitude, but of opposite signs. 

Using the equations satisfied by these fields we also obtain that: 

J A^3r]-'^) = J[iujei(^ A^3r'^)-iim{v^A*3V''% (52) 

D D 

and 

|d(^A*3r7+'°) = y"[^^eo(i+^A*3^+'°)-iuW(^A*3^+'0)]. (53) 

Adding, we obtain that 

j [iue(^ A *30 - A *3i^)] = j ^+'° A *3^+'°. (54) 

Br bBR 

The real part of the left hand side is non-positive as Re(ia;e) and — Ile{iZJJl) 

are non-positive and both ^ A *3^ and rj A *3r/ are point-wise non-negative. 
On the other hand we see that the radiation condition implies that the 
integral over the sphere can be written as 



^+.0 A ^3^7+'° = y ^ / A *3r?+'° + o(l). (55) 

bBR bBR 

As R tends to infinity, the real part of the right hand side tends to a non- 
negative number. 
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Arguing as in Miiller, we easily conclude that (^'^''^, ri^'^) vanishes identi- 
cally in ri. Equation ()49p then implies that j\ and mi are also zero. Finally, 
equations (jH]) , (jl2|) , and the orthogonality of the Hodge decomposition im- 
ply that the Debye sources, (rg, qo,ri, qi, j H^mn)-, are themselves zero. □ 

Now suppose that D = DiU - ■ ■UD]\i, has > 1 connected components. 
In the component, Di, we have EM parameters (eo«)^o/)- Oii each boundary 
component we have (inner and outer) currents {jii,'mii,jQi,mQi); on F; they 
satisfy the relations: 

Joi = \ —^dii and moi = . —Uimu l = l,...,N, (56) 
V <=o/ V ^oi 

where Ui : C^{bDi; A^) — )• C^{bDi; A^) is a choice of clutching map satisfying 
the conditions in Definition[TJ These boundary currents are therefore defined 
by Debye sources 

1^1 = (.roi,qoi,rii,qii,jHbiTT'Hl), l = l,...,N, (57) 
via relations analogous to ()4ip and ()42p : 



V 

mil = ^uJ{^/JIle^d^Roqll - m\ — *2 drRoqoi) + mm- 

V ^1 



(58) 



Now that we have introduced the relations between the Debye sources, and 
boundary currents, it is useful to modify the notation introduced in (p2|) so 
that the fields depend explicitly on this data. We use Vi^i to denote the 
currents {jn^mn) defined on F/ by ([58]) . and Vi^ the currents, {johi^oi) 
defined on this surface by (f58]l and (|56l) . For a; 7^ 0, we can identify 



J"=^(ei,/xi,a;,F/,Pi,i) = T^{ei, fii,u),Ti,jii,mii) and 
-^"^ (eo« , /^o« , F, , Pi_o) = -^"^ (eo/ ,m,^,'^i,joh '^oi ) ■ 



The fields F^{ei, jii^uj, F^, P^^i), and -F^(eoh /io«, F^, P^^q) are smooth func- 
tions of the Debye sources, even as uj goes to zero. If we need to refer to 
specific field components, then we use the notation F^,F^. 

The boundary condition, ()30p . defining the dielectric problem is now 
assumed to hold on each boundary component. We define a solution to the 
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problem via the following prescription: 

N 

{^,V) = Vj^+(ei,^i,w,r/,P/,i) in Q 

fit (60) 

(^,r/) = T'{eoi,noi,uj,Ti,'Difi) in A / = 1, . . . ,iV. 

The currents in these expressions are assumed to be defined in terms of the 
scalar sources and harmonic 1-forms by the relations in (j58p and (|56p . 
We can now prove the uniqueness theorem in this case. 

Theorem 2. Let bD have N components and fix u; ^ 0, so that 

{ujiii,uj€i,ujfj,oi,ujeoi : l = l,...,N} 

have non-negative imaginary parts and the ratios Re(/xoi/eo/) > 0. Let Di 
denote the Dehye source data defining solutions to THME(kQi) in Di, resp. 
THME(ki) in $7 as specified in (j60|) . where the boundary currents satisfy ([56]) 
and (j58|) . If{$,r]) satisfies ([30|) with j"^ = m™ = 0, then the Debye sources, 
{Vi} are also zero. 

Proof. As before, Miiller's uniqueness theorem shows that rj) = (0, 0) 
throughout M^. To show that the data itself is zero, we simply apply the 
argument in the proof of Theorem [T] to one component of at a time. Fix 
an 1 < Iq < N, and define 

ili„Vj=T+iei„^ii,,u;,r,Vi,^o) in \ Di, 

N 

(61) 

These fields satisfy THME(A;;o) in \ Di^ and THME(/ci) in Di^. Since 
(^, ij) vanish identically, by using the jump conditions and (j56|) . we see that, 
along Tig, we have the boundary data: 

?/ot = \//^*2 mii,^, {in.%)t = -V^Jiio, (62) 

and 

^iot = -\^*2l^rniig, {inVlg)t = V^^Jiio- (63) 
By applying the integration by parts argument from the proof of Theorem [H 
we can conclude that (^ig) ^i|^) vanishes identically. The boundary condition 
in (|63p then implies that 

{jii,,muo) = {0,0). (64) 
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Repeating this for each boundary component, and applying the Hodge de- 
composition completes the proof of the theorem. □ 

The naive limit of the dielectric problem, as w — t- 0, leads to a pair of 
uncoupled, underdetermined problems: 

= d*^ = and dr] = d*r] = 0, (65) 

where t?) satisfy (j30|) along bD. If we let denote a function which is 
harmonic in M'^ \ bD such that 

\bD= U~ \bD, (66) 

then the 1-form defined by 

= du^ in VL and = du^ in D, (67) 

satisfies the system of equations ()65p and the tangential component has no 
jump across bD. As we shall see, our representation of solutions suggests 
that, at cij = 0, we should also make use of conditions on the normal com- 
ponents of the form: 

eoC - ei^n = / and ^o(*3^?~)n - Aii(*3»?^)n = h. (68) 

If 7^ 0, then a condition of this type is a consequence of (|30p and the 
Maxwell equations. In the case of a scattering problem, 

/ = eijlf and h = ^imjf , (69) 

which can be applied for any w. If we append this condition, then the 
uniqueness theorem above extends to the zero frequency case. We state the 
result for the ^-field, the analogue for the r/-field follows by application of 

Theorem 3. Let ^ be an outgoing (zero-frequency) solution to (I65p defined 
in M'^ \ bD such that, along bDi, I = 1, . . . , N, we have: 

$t - = and eoiC - ei^+ = 0. (70) 

Suppose that {ei,eoi '■ 1 = 1,..., N} all have positive real part. If the normal 
components have mean zero on every component ofbD, then the solution 
^ vanishes identically. 
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Proof. For a a closed 1-form along bD, we let [ctJdR denote its class in 
H^^{bD). As d^^ = 0, the tangential restrictions are closed 1-forms. 
The tangential boundary condition in (j70p implies that 



[itU = [^t"]dR. (71) 

On the other hand, the Mayer- Vietoris theorem implies that [^t^]dR and 
[^r]dR belong to complementary subspaces of H^^{bD) and therefore both 
must vanish. Hence there are functions defined on bD so that = df^. 
We let u~ denote the harmonic function in D with u~ \bD= /~; and 
the outgoing harmonic function in Q with n"*" \bD= f^- The field ^~ — du~ 
has vanishing tangential components on bD and is therefore identically zero 
in D. The tangential components of — du^ are zero, showing that this 
difference is a sum of Dirichlet fields. As has vanishing mean on each 
component of bD it follows from Theorem 5.7 in [2] that, in fact, = du^ . 
If we let VLh = Oni?/j(0), then the normal boundary condition in (|70p easily 
implies that 

TV 

ei j \du^\^dx + ^eQi j \du-\^dx= Ju+d^dS. (72) 

Qn '=1 D, bBn 

As ti+ is outgoing, the limit, as i? — )• oo, of the integral over bBji is zero. 
This completes the proof that ^ = 0. □ 

The uniqueness of the representation of solutions to (|70p at zero fre- 
quency, via Debye source data turns out to depend on the choice of clutching 
map U : 'H^{bQ) — )• 7^^(60). Of course U = -^2 is such a map. Whether this 
suffices to prove uniqueness at cj = turns out to depend on a surprisingly 
subtle property of the Hodge star-operator on bD relative to the splitting of 
H^^{bD) into the disjoint Lagrangian subspaces (w.r.t. the symplectic form 
defined below in ([76]) ) 

Hl^{bD) ^ Hl^{D) \bD ®Hl^{n) kD . 

We let 'H\){bD) denote harmonic representatives for the image of the injec- 
tive map H\^{D) ^ H\^{bD), and 7iQ{bD) harmonic representatives for 
the image of H^j^{Q,) H^^{bD). As noted above, the Hodge star-operator 
maps harmonic forms to harmonic forms. 

Definition 2. We say that bD is H-generic if 

*2nh{bD)nnUbD) = {o}. (73) 
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Evidently a simply connected manifold is H-generic. More generally, 
the property of H-genericity depends only the conformal structure on bD 
(induced from its embedding into M^) and the splitting 

Hl^ibD) = Hl^{D) \,D ®Hl^m ItD, (74) 

which depends on the isotopy class of the embedding of bD ^ R^. 

It is a deep theorem in algebraic geometry, stemming form work of Rie- 
mann, that the set of H-generic structures is the complement of a real ana- 
lytic hypersurface in Teichmiiller space, and therefore open and dense. On 
the other hand, there certainly exist surfaces for which (j73p fails. For exam- 
ple if bD is a torus of revolution, as described in Example 1.5 of [3], then we 
have -k2'H\i{bD) n 'H]j{bD) = 'H]-,{bD). If we were to take U = -k2, in such a 
case, then the Debye source representation at zero frequency would have a 
non-trivial null-space of dimension equal to dim*2^^(6D) r\'}i\){bD). This 
would then lead to a mild form of low frequency breakdown. To avoid this 
eventuality we need to make a different choice of lA. 

Definition 3. A clutching map U : ^{^{bQ) — )• ^{^{bQ,), is admissible if, for 
each component bDi of 60, we have 

Ui *2 nhXbDi) n *2nUbDi) = {o}. (75) 

It is easy to see that such maps always exist. As we do the construction 
one component at a time, we can restrict attention to the case that bD is 
connected. We let uj denote the hermitian symplectic form 

uj{a,l3) = y a A/3. (76) 

bD 

This is a non-degenerate pairing that is well defined on H^^{bD). Moreover, 

uj{*2a, i.2P) = uj{ct,P). (77) 

The topological interpretation of this pairing implies that both *2^£i(6-D) 
and *2'HDc{bD) are Lagrangian subspaces; that is, uj restricted to these 
subspaces is identically zero. In fact 'H\){bD) is a Lagrangian subspace 
complementary to 'k2T~i\){bD) , and these subspaces are complexifications of 
real Lagrangian subspaces. The map lAi = Id on 'H^{bDi) is always an 
admissible clutching map. For if 

aG*2'Hh,ibDi)n*2'Hhf{bDi), (78) 

then 

*2cxenUbDi)nn'Dc{bDi) = {o}. (79) 
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Theorem 4. Suppose that every component of bD is H-generic, in which 
case we takeU = or thatlA is an admissible clutching map. If an outgoing 
(zero-frequency) solution to (j65p in \ bD, defined by Debye source data 
via (j60p . vanishes, then the data defining it vanishes as well. 

Proof. Note first that, at zero frequency, the Debye source data uncouples, 
with ^ determined by 

{{rounumm) : l = l,...,N}c Mr,o x nHbD) 

via (f60]) . and ij determined by 

{{qoi,qiiJm) ■■ l = h...,N} cMr,oxnHbD). 

We restrict our attention here to the electric field and suppose that the 
Debye source data determines a 1-form ^ that vanishes identically. First we 
assume that the harmonic 1-forms {rriHi : 1 = 1,... ,N} are all zero. This 
is always the case if bD is simply connected. In Di we have 

r = -^/mdct>T, (80) 

where 

9Q{x,y)roidA. (81) 



bD 



The fact that dcl)^ = implies that (f)^ is constant, which implies that <f)'^ , 
its continuation to Df, is also constant along bDi. Integrate by parts, using 
that dn(j)~^ = rQi, and 0/ = c a constant, to obtain 

\V^l\^dx = c j roi. (82) 

Df bDi 

Either because c = or because r^i has mean zero we see that, in fact, 
0^ = as well, which implies that tq/ = 0. Similarly, in we see that the 
harmonic function 

N 



j 9o{x,y)riidA (83) 



is constant and therefore 0. This function is continuous across each boundary 
component, so it is also zero in each component of Di. The jump in the 
normal derivative of 0/ across bDi is rn, which must also vanish. This 
completes the case where mn = 0. 
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It remains only to show that the harmonic 1-forms {rriHi} vanish. For 
each 1 < Iq < N we use the prescription in (j6ip to define a harmonic field 
in \ bDig . The relations in (|62p and (j63p show that along bDi^ 

'^lot = \//^*2 rnnio and = -^-k2 ^rnnio, (84) 

where mHi^ is a harmonic 1-form on . As d^^^ = in M'^ \ hDi^ , we see 
that 

*2 G -^^k (6Ao) and -k2UmHif, G ?^l)c (^Aq)- (85) 

'o 

The assumption that lA is an admissible clutching map implies that *2i^Hio — 
0, completing the proof. □ 

To conclude this discussion, we show that in the case that bD is not 
H-generic, and we take U = -k2, then the Debye source representation has a 
non-trivial nullspace. Suppose that bD is connected and 

*2 mn G ^nlibD) n V.'oibD). (86) 

We let j = 0, 1, denote the fields defined in \ bD by 

= J-±(ei,/ii,0,6D,{0,0,mH}) ^ ^ 

, ^ (87) 
^0 = -^f (eo,Aio,0,6-D,{0,0,*2m//}). 

The assumption in (|86l) implies that ■^2i^H G 'HjjibD), and m// G ^^{bD). 
Arguing as in the proof of Proposition 7.12 in [4j, we see that the cohomology 
classes [$,1] and [^q ] are trivial. Thus we can find unique scalar Debye 
sources (ro,ri) defining fields dcf)^ so that 

d(l)'l + = and d</>o + = 0. (88) 

This shows that the null-space of the Debye representation at zero frequency 
includes a subspace isomorphic to -k2T~i^{bD) n 'H]-,{bD). From the proof 
of Theorem H] we easily conclude that this is exactly the null-space of the 
representation. 

4 Integral Equations for the Dielectric Problem 

In this section we derive Fredholm integral equations of the second kind for 
solving the dielectric problem, in terms of the Debye sources. The dielectric 
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interface boundary conditions are formulated in ()30p . We rewrite them 
slightly for this section 



*2 i^i'' - = MjT and *2 [(*3r7+'^)i - (^ar?"'")*] = ^m}^. (89) 

The data (j™,m™) are arbitrary 1-forms defined on F. We usually assume 
that these are defined as boundary data of a solution (^™,'/7™) to the time 
harmonic Maxwell equations m.Vt. To obtain equations in terms of the Debye 
sources, we apply the operators Go *2 and dr*2 to these equations. In 
the simply connected case this suffices to solve the original problem, as a 
1-form fc on r is specified by the scalar functions •k2drk and d^k. 

The tangential, and normal operators, 7^ (k), (k), and (k) , ( /c) , 
used to construct ^^'^ and rj^'^ are defined in (|23p and (j24p . There are differ- 
ences between these operators and those used in the perfect conductor case. 
For the latter case, we used (jlSp to define j in terms of of a single set of De- 
bye sources. Here, for Z = 0, or 1, the vector sources J;,m/ in (f23]l and (i24l) 
depend on four scalar potentials {rQ,qo,ri,qi), and possibly the harmonic 
1-forms jfj,mH, which we abbreviate as {r,q), {r,q, j jj^uih) resp. More- 
over, the vector sources jf\ and mi are independent of one another, with the 
inner and outer sources, (j'oj'tio)) satisfying the relations in ([Mil- 

Notice also that the wave numbers are different in different regions of space, 
so that, e.g. we use 7^~(A;o) to define the ^-field in the interior of D, and 
T^{ki) to define the ^-field in the exterior domain il. Finally, it should be 
noted that, when hD has several components, the dependencies between cur- 
rents and scalar Debye sources stated in (i58]l are local within each boundary 
component. This is the case even though the formulae in (j60p imply that 
the integral equations we ultimately have to solve intertwine the different 
boundary components, albeit via smoothing terms. 

While our approach was motivated by that of Miiller [7], it is worth 
noting that there are two important differences. He used two vector fields 
(jfi,mi) as unknowns combined with a relation of the form 

Jo = \ —Ji and mo = W — mi 
V eo V /"o 

and proved well-posedness for the resulting integral equation. In order to 
avoid low- frequency breakdown and achieve graceful uncoupling of the fields, 
the introduction of the Debye source representation necessitates a more del- 
icate analysis of the interior and exterior representations, and the use, in 
some cases, of a non-trivial clutching map. 



21 



4.1 The Single Component Case 



We first write the integral equations for the case that bD is connected. In 
this case {ro,qo,ri,qi) G -A4p g ^^"^ {3h^''^h) G %^{bD) x l-0'{bD) need 
to be determined. When possible, we omit the explicit mention of these 
arguments below. Our first two equations are 



Go4[(^^ -^7' )] = G^d'^r and Go4[(*3^+'')t - (*3^"'°)t] = G^d'^m}-. 

(90) 

More explicitly, these equations are 



Go *2 d^[^/e^Tr|~ (h) - y/e^Tr, {ko)] = God^m' 
We use the fact that 



(91) 



GodfdrGfc = ^+0(-l), (92) 



to deduce that these equations are of the form: 



ri + ^—ro + 1. o. t. = Go^rJ 



4 4 



qi - -^go + 1- o- 1- = Goctr"^ 



(93) 



Here and in the sequel, 1. o. t. refers to pseudodifferential operators of nega- 
tive order (when applied to {r,q)). 

To derive the other equations we observe that for a 1-form a smooth up 
to r we have that 

driat) = {da)t, (94) 

and 

d$^'^ = iujfj.ir]^''- and r?^'' = -iuei ^3 (95) 
This means that 

(d^±'')t = iojfii{r]^^% and {d *3 r/±'')t = -i^^ez(*3^±'')t . (96) 

The right hand sides of these equations are essentially the normal compo- 
nents for 77^'' and respectively. On the other hand, (|94|) implies that 

idt'')t - {dr% = drr and (d*3 r/+'^)t - (d*3 = drm'^, (97) 
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and therefore it is reasonable to assume that there are functions /, h defined 
on r, of size 0(1), so that 



drj" = iujfdA and drm}'' = iujhdA. (98) 

This is automatic if the data (j™, m™) come from a solution of the THME(A;i). 
Note that (/, /i) E A^r,o- 

With this understood, the second set of equations becomes: 



-eiy/JIiM^{ki) + eQy/JI^M^ {ko) = h. 
These equations are of the form: 



(99) 







2 + 


2 




\/W)eo 


2 + 


2 



qoj +l.o.t. = /, 
ro I + 1. o. t. = h. 



(100) 



These relations, along with (I93p show that (I9ip and (I99p are a Predholm 
system of the second kind for {ro,qo,ri,qi), provided that 

MoMieoei / and (^uq + m)(eo + ei) / 0. (101) 

If the genus of F is p 7^ 0, then we also need to determine {jjj^rriH) G 
?^^(r) X V}{T). Equations ^ and ([Ml) do not quite suffice as the right 
hand sides of these equations explicitly annihilate the harmonic projections 
of J™ and m™. We have two distinct choices as to how we should augment 
these equations to capture the projection into the harmonic vector fields. 
On the one hand we could use a basis ^ = {tpi, . . . ,'ip2g) for ^^"'^(r), the 
harmonic one forms on F, and then augment the equations (j9ip and (|99p 
with the 4(7 equations: 

([(*3r7+'^)t-(*3r7-'%)],^) = (m-,M/). 

We can also follow a more traditional, and geometric approach. Let 
{{Aj,Bj) : j = 1, . . . , 51} be a basis for Hi{T), normalized as above: the 
j4-cycles span [H^{D)]^, and the i?-cycles span [H'^[Q)]^. We can then 
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augment the equations ([9T]) and ([99|) with conditions on the circulations: 



i^t'' - ^?'] = J r 

j [(*3r7+'')t - (*3r7-'°)0] = / 

■Aj Aj 

j [(*3r/+'^)* - (^sr?-'")*)] = / m" 



(103) 



If we have Debye source data in the null-space of ([9T]) and ([99|) . then the 
tangential fields [^^'^ — ^7'°] [(^s^?''''"^)* — (^s*?"'*^)*)] are harmonic 1- 
forms. As H\^{bD) and Hi{bD) are dual vector spaces, via this pairing, 
such data lies in the null-space of (jl02p if and only if it is in the null-space 

of (umi). 



4.2 The Multiple Component Case 

In this subsection, we assume that bD = bDi U • • • U bD^^f. Taking account 
of the differences between (I33p and (I60p . we see that the equations in this 
case differ somewhat from the equations in the case where bD is connected. 
The differences are all in the form of smoothing operators, and therefore 
the computations of the leading order terms apply mutatis mutandis. For 
clarity, we augment the notation for the tangential and normal restric- 
tions to specify a particular component of the boundary, e.g. it should 
be understood that {ki,bDi) is a tangent field on bDi that depends 

on (,roi,ru,qQi,qii,jfji,mHi). For m / /, we use J"+/(ei, /^i, w, F^, 

J^^'J^{ei, fii,uj,rm,T^m,i), to denote the tangential and normal components, 
resp., of the ^-field along bDi of the fields defined by the sources (j'lm) ''^im) 
on bDm- The notation F^'l,J-^^n has the analogous meaning for the rj-field. 

With these notational conventions, we can now give the equations for the 
multiple component case. We let (jf™, m™) denote the incoming fields along 
bDi. As in the previous case we assume that these fields arise as restrictions 
of a single solution (^™, ?7™) of THME(/ci), and therefore, there are functions 
{fi,hi) of size 0(1) so that, along bDi we have: 

dY^Jf = iujfidA and dr^mp = iujhidA. (104) 
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Now fix 1 < / < A^. Along bDi the tangential equations for the solution 
of the dielectric problem in terms of the Debye sources reads: 



Go *2 dri^^Trj' {h) - ^/e^iTr, (%)] + 

God^ ^ ^ J^^^t (ei , /ii , u},Tm, 

Away from F; we have the relations 

dT^{ei,fii,Ti,'Di^i) = itjfiiT:^{ei,fii, 11,7)1^1) 
d-ks F:^{ei,^i,Ti,Vi^i) = -iujei +3 J"^(ei, /ii, T;, 

These relations and pU4|) show that the normal equations are 

fIl^/el^^^{kl,bDl) - iiQi^/e^iNr,{kQubDi)+ 

/^l -^T?in('^li /^li'^irm, Vm,l) 



(105) 



(106) 



fl 



*2drf 



if w ^ 



-eiy/JI^Af^{ki,bDi) + eoiy/JmM^ {koi,bDi)+ 



= ,, = (5*Z!^,„^„). 

(107) 

As noted above, the new terms in these equations are smoothing operators. 
Thus, with the obvious changes in notation, the equations for the leading 
order parts, ([93]) and (jlOOp apply equally well in the multi-component case 
to show that these are again Fredholm equations of the second kind for 
{{i^ohi^ih QohQii) ■ I = I, ■ ■ ■ , N}. To complete the system we append either 
the 4p equations in p02p . or those in p03|) . where p is the total genus of F. 
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4.3 Low Frequency Behavior 

Suppose that we are given Debye data in the null-space of the system defined 
by (fTOSD . (fTOT]) . and (fT02]) . If a; / 0, then it is apparent from the Hodge 
theorem that {^^,rj^) satisfy ([5U|) with = = 0. Hence Theorem [2] 
shows that the Debye data is in fact zero. If we are given Debye data, so 
that, at cj = 0, the limiting homogeneous boundary conditions ([70]) hold, 
then, provided that U is an admissible clutching map, Theorem [J] implies 
that again, the solution (^^,'?7^) is identically zero, as is the Debye data. 
That the niecin valuGS of the nornicil components i^in^'^f^n 

*3 r/) vanish on 

each component of bD follows immediately, as this is true for all solutions 
defined by Debye data. 

Thus, if we make the assumption that the incoming field satisfies the 
conditions in ()104p . and U is an admissible clutching map, then this system 
of equations displays no low frequency breakdown. The uniqueness results 
for the Debye representations. Theorems [2] and show that, using the ansatz 
in (j60p . the only Debye data which leads to a solution to the homogeneous 
dielectric problem is zero. As (llOSp . ()107p . and (I102p is a system of Fredholm 
equations of second kind, this proves the solvability. 

Theorem 5. Let D be a union of smooth bounded regions in M'^ with con- 
nected complement Q. Suppose further that U is an admissible clutching 
map. 

1. //w / and we assume that 

{uJHi,ujei,ujfiol,ujeoi : l = l,...,N} 

have non-negative imaginary parts, and the ratios Re(^o«/£oz) > 0, 
then the system of equations (jlOSp . ()107p . and ()102p for 

{{roi,rii,qoi,qu,jHh'n^Hl) ■ l = '^,---,N}, 
is solvable for any right hand side (j™,m™). 

2. If uj = 0, then, assuming that, for some (/), {e^'^ei^e^'^e^i : 1 = 1,..., N} 
all have positive real part, the ^-equations in (I105j) . (jl07p . and (]102p 
for the Debye sources {{rQi,rii,mj{i) : 1 = 1,... ,N}, are solvable for 
arbitrary admissible data dfj"^, {/i)---)/Af} ".n-d Here data 
is admissible if drj"^ = and 

J fidA = 0. (108) 
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3. The analogous result holds for the r]- equations assuming that there is 
a 9 so that {e*^/ii, e'^/ig/ : / = !,..., all have positive real parts. 

Remark 2 . The positivity conditions in this theorem are true if Im a; > and 
the physical constants satisfy the conditions given in Q. As noted above, 
for an H-generic boundary we can let U = •k2- H-genericity is an open and 
dense condition on surfaces in M^. 



5 Low frequency behavior of the perfect conduc- 
tor problem 

The operators defining the Hodge star (*2) of the tangential components of 
(^^,*3r7^) are given by 



Tfik) = —--Kir + ikK2,t j - K^m 



2 

T^{k) = ^-Kiq + zkK2,t m + K^j 
The equations for the normal components of the rj-fields are 



(109) 



^^{k) = -K^q + Ksj + ikn, ■ Gk m (110) 

For the perfect conductor case, we introduced a single pair {ri,qi) of 
scalar Debye sources for each component, Ti, of bD and a single harmonic 
1-form jfji- The currents {ji,mi) are given by 

= ik{dTiRori - -k2dTiRoqi) + 3 m and mi = *2jV (HI) 

Using {T^ {k){r, q, j jj) , {k){r, q, j fj)) to denote these operators restricted 
to A^(r) X H^(T), with these relations, the hybrid system of integral oper- 
ators is then defined to be 




The range of *2drT^{k) is contained in the space of functions on T with 
mean zero on every component. 

As noted in for cj / 0, the nullspace of this system consists of data 
such that is a harmonic 1-form. We can therefore append relations of 
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the type given in ()102p or (jl03p to get an invertible system. A system of 
the type in ()102p . however, has serious conditioning problems as w — >• 0. 

If we suppose that the incoming data (^^,^7™) is a solution of the 
THME(A;), then an auxihary equation similar to that in (jlOSp can be em- 
ployed, which does not suffer from conditioning problems at k = 0. We let 
{Aj, Bj : j = 1, . . . , g} he a basis for Hi{T), normalized, as above, so that 
the {Aj} are a basis for [Hl^{D)]^ and the {Bj}, a basis for [//jf^(O)]^. We 
augment 



r 




(113) 



with the equations 



*2r+(fc) = J c (114) 
i/*.r/w = i/c. (115) 

Bj Bj 

Theorem 6. For a; / 0, the null-space of the full system (|113p . (jll4p . 
and piSp is trivial. 

Proof. If {r,q,jj^) is in the null-space of (|113p . then the tangential 1-form 
^j' £ n^iV). If the integrals in (fTTil) and (fTT5|) vanish, then is the 
harmonic representative of the trivial cohomology class in H^^{T), and is 
therefore zero. Hence = 0; as uj ^ 0, Theorem 7.1 in f3] shows that the 
data {r,q,jjj) must also vanish. □ 

At w = the current j = Jhi we see that 

if = ^ - dvGQr + ★2i^4 *2 3h 
[*3»?^]t = 1^ - drGoq - -k2Ki3H- 

Since the cycles {Bj} span [//j^j^(r2)]-'-, there are smooth surfaces {Sj} con- 
tained in VL so that bSj = Bj. Thus, using Stokes theorem and the equations 
di'^ = ikrj'^, d^™ = ikr]™, we can rewrite the conditions in (jllSp as 

j M+{k){r,q,3H) = j rf-. (117) 
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The advantage of this formulation is that it has an obvious smooth hmit as 
A; — >• 0. In fact, replacing (jll5p with these equivalent conditions allows us to 
verify that the augmented system of equations (|113p , (|114p , and (jll7p does 
not suffer from low frequency breakdown. 



Theorem 7. // {r,q,jfj) G Air x 'H^(r) defines solutions {^^,ri^) to 
THME(O) such that 



and the integrals on the left hand sides in ()114p and (I117p vanish, then 
r = q = 3h = 

Proof. If, for /c = 0, ^"'^ comes from data {r,q,jff) G Mt,o x ^^(T) that is 
in the null-space of ()113p . then, also using the fact that d^^ = 0, we see 
that satisfies 



Thus is a harmonic 1-form that belongs to the image of H^^{Q) 
H^^{bQ). If the integrals in the first line of (jll4p are also 0, then this implies 
that the cohomology class represented by is trivial and therefore = 0. 
The integral of in$^ vanishes over each component of T and therefore the 
standard uniqueness result, see p], for outgoing harmonic 1-forms in shows 
that $,+ = 0. 

If T]^ is in the null-space of (jll3p . then the normal component, rj^ \bn= 
0. This implies that rj^ is a harmonic Neumann field, which therefore is 
determined by its class in H^-^{Q,,bQ). As the surfaces {Sj} define a basis 
for H2{^l,bCl) = [H'^^{Q,,bQ,)]' , this class is, in turn, specified by the values 
of the integrals in (I117p . Hence if these integrals are also zero, then t]^ = 
in Q as well. 

This shows that if there is data {r,q,j}{) for which both 



and the integrals in the auxiliary conditions, (|114p and ()117p . vanish, then 
(^^,■77+) = (0,0). From the jump relations implicit in (|116p we conclude 
that 




(118) 



(119) 




(120) 



Jh and [★3/7 ]t 



(121) 
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Since d$,~ = d-k^ r]~ =0 this implies that both jfj and *2J_h' must belong 
to Im(ffjjj^(D) ^ H^^{bQ,)). This is a Lagrangian subspace with respect to 
the wedge-product pairing, and therefore 

J//A*2jH = 0. (122) 



bn 

Hence jjj = as well. The fact that (^^,^7^) = (0,0) now implies that 
their normal components vanish. Coupled with the vanishing of jjj we see 
that this shows: 

^ + i^o(O)) q = + i^o(O)) r = 0. (123) 

It is a classical result, see [21 [H], that these operators are invertible and 
therefore (r, q) = (0, 0). This completes the proof of the theorem and demon- 
strates that the hybrid equations along with the auxiliary conditions in (|114|) 
and (|117p do not suffer from low frequency breakdown. □ 

For any harmonic 1-form jfj G 'H^{r), we see that, at w = 0, 



^ + *2^4 *2 Jh e Im(/fiR(0) ^ Hl^iV)) 

^ - *2KdH e im(//dR(^^) ^iR(r)). 



(124) 



So these inclusion maps have rank equal to ^ dimff^j^(r). The system (IllSp . 
(jll4p . and (jll7p splits into two almost uncoupled systems, one for as a 
function of d^^f and the integrals in (I114p . and the other a system for r]^ 
in terms of ij and the integrals in ()117p . These systems are coupled only 
through jjj, which appears in both rows of (]113p . 

The coupling can be effectively removed by choosing a basis {ipj : j = 
1, . . . , 2g} for H^ir), so that {ipi, . . . , ipg} spans the nullspace of the operator 
*2j//— 2*2-f^4j_f/) ^■s a ™ap from ^{^{T) to H^j^iT), and {ipi+g, • • • , ip2g} spans 
the nullspace of {jj^ -|- 2 *2 K4 *2 Jh) • ^^(r) ^dRi^)- "^^^ uniqueness 
theorem just proved, and the invertibility of r ^ G^r on C°°(r), show that 
these subspaces are complementary, and therefore 

y + *2i^4*2 j : i = 1, • • • , 5} (125) 
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spans the image of -ffdR(^) ill ^dR(r)- Using the conditions in (fTTil) we can 
obtain coefficients {aj} so that 



^ = + E "i ( T + ) V', (126) 

is trivial in H^-^{r). The Debye source r is then found by solving the equa- 
tion: 

- GodrdrGor = God^K. (127) 

An analogous discussion applies to the equations for the magnetic field. 

For numerical purposes it is possible to split the difference in the inte- 
grals over the S-cycles. The integrals of the data 77™ are computed over 
the spanning surfaces {Sj}. The integrals of over cycles {Bj} can be 
accurately computed to order k, by taking advantage of the fact that when 
A; = 0, we have 

J *2r+(0) = J dt=0 for j = l,...,g. (128) 

Bj Sj 

We can therefore rewrite the equations in ()115p as 

1 1 *2[r+(fc)-7^+(0)]=i J r/-. (129) 

The accuracy of the calculation on the left hand side is retained by observing 
that, as \k\ is assumed to be small and \x—y\ is bounded on T, we can employ: 

^[..(x,y)-.o(x,.)]=^El^^^ (130) 

to avoid having to explicitly divide by k, which thereby avoids catastrophic 
cancellation. 

To use either (jll7p . or (|129p . requires finding surfaces {Sj} that span 
the 5-cycles, i.e. 

bSj=Bj for j = l,...,g, (131) 

and performing certain 1- and 2-dimensional integrals. The use of (jl29p 
has the advantage that the integral over the (artificial) spanning surface Sj 
involves a known and typically smooth incoming field, while the complicated 
integrand [7^"^ (A;) — 7^'^(0)] needs only be computed over Bj which lies on 
the surface bD itself. 
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6 The Dielectric Problem on a Torus of Revolu- 
tion 



In the case where T can be described as a surface of revolution, we can 
represent the unknown charges ri, ro, qi, and qo in terms of their Fourier 
expansions in the azimuthal variable. More precisely, we assume that the 
surface F is given by 

x(t,9) = p{t) COS0 
y {t, 9) =p{t) sin 9 
z{t,9) = z{t), 

where the generating curve 7(t) = {p{t),z{t)) is smooth, with t G [0,L] and 
9 S [0, 27r] (Fig. [l]). It is straightforward to verify that 

/ f — sm9 cos 9 \ , . / / ^ . ^ z'(t)\ 

form a basis for the harmonic 1-forms on the surface F. 




Figure 1: A surface of revolution F (of genus 1), defined by a generating 
curve 7(t) = {p{t),z{t)) which is rotated about the z-axis. We let p{t) = 
2 + (1 + 0.2 cos(4f)) cos(t), z{t) = 2 + (1 + 0.3 sin(4i)) sin(t). On the left is 
plotted the intersection of the surface with the xz-plane. On the right, the 
the z-component of the harmonic 1-form j ^2 is plotted as a color value on 
the surface. 

In this setting, we may expand each scalar Debye source ri and qi, for 
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/ = 0, 1, as a Fourier series in the 6 direction: 



ri{t,e)= J2 n,n{ty''' 

n=—oo 
oo 

qi{t,9)= QiAty'"'- 

n=—oo 

Using this representation for the scalar sources leads to our electric and 
magnetic current-like variables via equations ()4ip and (j42p : 



'^^""^oo (132) 



where r and are globally defined orthonormal unit co- vectors defined along 
r. Recall from above that the harmonic 1-forms j and j ("^Hi and 
mH2) are constant in the azimuthal direction, so that they couple only to 
the purely axisymmetric (n = 0) mode. 

In short, the linear system ([9T]) . ([99]) . (jl03p can be solved for each az- 
imuthal mode separately. For nonzero azimuthal modes, where the harmonic 
1-forms play no role, only equations (|9ip. (j99p are needed to determine 
i^i,n{'t),Qi,nit) for / = 0, 1. For each mode the set of equations corresponds 
to a Fredholm system of equations of the second kind along the generating 
curve. 

In a companion paper [5] , we have developed a high-order accurate solver 
for scattering from closed surfaces of revolution, with a detailed description 
of the full algorithm. Here, we simply note that the method of [3] assumes 
an equi-spaced discretization of 7(t), and uses a pseudospectral approach 
for applying (or inverting) the surface Laplacian, the surface gradient, and 
the surface divergence. That is to say, differentiation is carried out in the 
transform domain, while multiplication by a variable coefficient function is 
carried out pointwise along 7(t). The principal value and weakly singular 
integrals which appear are computed using generalized Gaussian quadrature 
rules [Tj with eighth or sixteenth order accuracy. 

One issue that requires some care is the mean zero condition on the 
scalar sources. We showed, above, that the integral equation is invertible 
as a map from mean zero functions to mean zero functions. This condition 
is automatically satisfied for non-zero azimuthal modes (which necessarily 
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integrate to zero). For the zero mode, one can either solve the integral 
equation iteratively, relying on the fact that the operator projects onto mean 
zero functions or add a rank-one modification to the system matrix to enforce 
invertibility of the discretized system and solve it directly. We have chosen 
the latter approach, described in more detail in [S]. 

Remark 3. In the case of a perfect conductor, there are only two scalar source 
functions {r{t, 6), q{t, 6)), expanded as above, with the electric and magnetic 
current-like variables computed from (jlSp and the condition m = *2 J • 

Remark 4. In the dielectric case, we still need to specify the choice of a 
clutching map. Our analysis above showed that a sufficient condition for 
uniqueness that avoids low frequency breakdown is to set U = Idon T-i^ihD). 
Fig. [2] displays the condition number of the linear system for various choices 
of Lo and lA. 

Remark 5. The second kind Fredholm system derived earlier for scattering 
from a dielectric is free from low frequency breakdown, as well as spurious 
resonances. An illustration of this is contained in Fig. [3l which shows the 
relative P error in the electric and magnetic fields for the n = mode in the 
exterior and interior of the object in Fig. [TJ The surface was discretized at 
200 points in the t direction using a sixteenth order hybrid Gauss-trapezoidal 
quadrature rule (see [1]). In order to calculate the accuracy of the scheme, 
two extra surfaces of revolution were constructed - one inside the scatterer 
and one outside. Smooth generalized Debye sources and harmonic vector 
fields were specified on each, and the electric and magnetic fields generated 
with the appropriate material parameters were evaluated on the surface 
of the scatterer. The source surface inside the scatterer gives rise to a 
field which is valid in the exterior of the scatterer, and the source surface 
outside the scatterer gives rise to a field which is valid in the interior of 
the scatterer. The tangential and normal components of the difference of 
these fields were taken as the boundary data. Using this boundary data, 
the generalized Debye sources and harmonic vector fields on the scatterer 
were calculated using the integral equation formulation for the dielectric. 
Lastly, the scattered field was evaluated and compared with the known fields 
generated by the two source surfaces. 

Note that the error becomes even smaller at very small frequencies - this 
is because discretization errors that are introduced in the integral equation 
in terms that are 0{uj) and 0{uj'^) are suppressed as cj — >• 0. Only the accu- 
racy for the n = mode is plotted. The relative error, as a function of the 
frequency cj, in the n = mode is indicative of that in a fully reconstructed 
scattering problem (i.e., after reassembling the orthogonal, azimuthally de- 
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coupled, field components). 



Condltlan number vs. clutching map Condhlori numbei vs. trequency 




D 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1e-16 1e-14 1e-12 1e-10 le-OS 1e-a6 0.0001 0.01 1 

t Omega 



Figure 2: The plots above show the conditioning of the linear system for 
the purely axisymmetric mode depending on frequency and the choice of 
clutching map for the surface shown in Fig. [TJ On T-i^ibD) we define 
lA = cos t Id + sin t *2- On the left is plotted the condition num- 
ber for various values of t. Note that only when t is near 7r/2 and uj is near 
zero is the system ill conditioned. On the right is plotted the dependence of 
the condition number on oj when U = ^2- In both plots ei = 1.30, fii = 0.83, 
eo = 0.90, and /ii = 1.10. 



7 Summary 

We have developed a new integral representation for the solution of the time 
harmonic Maxwell equations in media with piecewise constant dielectric per- 
mittivity and magnetic permeability in M^. In the simply connected case, we 
rely on four scalar densities (charge-like variables) rather than surface vector 
fields as our unknowns, and obtain a system of Fredholm equation of the 
second kind. In the multiply connected case, we supplement these unknowns 
with a basis for the surface harmonic 1-forms. The principal advantage of 
our approach is that it avoids the low frequency breakdown inherent in the 
classical method, due to Miiller. Some subtlety arises in the selection of the 
interior and exterior representations. In Miiller's equation, the currents used 
to represent the interior and exterior fields are chosen to be scalar multiples 
of each other. In our case, we use a Hodge decomposition of the 1-forms 
and have introduced the clutching map in Section 3 to relate interior and 
exterior variables. We show that for stability in the zero frequency limit, 
the clutching map should act differently on the harmonic 1-forms and on 
their orthogonal complement. 
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Error in E and H 



- Interior E 
Interior H 

- Exterior E 
Exterior H 




4 5 6 

Frequency 



Figure 3: The plot above shows the accuracy obtained in the exterior and 
interior electric and magnetic fields at various frequencies for the purely 
axisymmetric mode. The material parameters were chosen to be ei = 1.30, 
fii = 0.83, eo = 0.90, and = 1.10. 



A disadvantage of our approach is that we must construct current-like 
surface 1-forms through a non-local operator. In the present paper, we 
do this by inverting the Laplace-Beltrami operator, but other (spectrally 
equivalent) procedures can be applied that do not require the solution of a 
linear system of equations [HI \5\ . 
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A The Boundary Operators 

We now describe the limits of the tangential components of ^ and 77 along 
r. Along r we can introduce an adapted local basis of orthonormal 1-forms, 
{ijOi,L02,i^}, and n the outward unit normal vector. In this basis we have 
= aj-coi + bj-UJ2 + c±v. The tangential part of at x € F is the 1-form 
$,^{x) restricted to directions tangent to F, i.e., T^T. In terms of components 
along F we identify the tangential part with 



If the 2-form, r] = e±uji A + f±uj2 A 1/ + g±oji A ljJ2, then the tangential 
components of 77^ along F are identified with the tangential components of 
the one-form i-nT]^ wedged with —v : 



1988. 



= a±uji + h±U2. 



(133) 



rjf = e±uji Av + f±0J2 A v. 



(134) 



The normal component of ^ is and that of r) is simply rj \y ■ 
The various boundary operators are given by: 




Kq is an operator of order —1 and Ki is an operator of order 0. 
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K2,n\j]{x) = in 
K2,t[j]{x) =-k2 



j^gk{x-y)j{y)-dxdA{y) , 
gk{x-y)3{y) -dxdAiy) 



These are operators of order —1. 

^3[j](a;) = d^gk{x - y) ■ {j{y) x n{x)) dA{y) , 

dgk 



dxQkix - y) {j{y) ■ n{x)) 



(x - y)j{y) ■ dx 



dA{y) . 



Ks is an operator of order and K4 is an operator of order —1. 
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